
volScalarField T
(
  IOobject
  (
     "T",
     runTime.timeName(),
     mesh,
     IOobject::MUST_READ,
     IOobject::AUTO_WRITE
  ),
  mesh
);

dimensionedScalar k 
(
  "k",
  dimensionSet(0, 2, -1, 0, 0, 0, 0),
  scalar(4e-8)
  // scalar(4e-5)
  // scalar(50)
);
